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ABSTRACT 

We extend our initial study of diffusive nuclear burning (DNB) for neutron stars (NSs) with Hydrogen atmo- 
spheres and an underlying layer of proton capturing nuclei. Our initial study showed that DNB can alter the 
photospheric abundance of Hydrogen on surprisingly short timescales (10^~'*yrs). Significant composition evo- 
lution impacts the radiated thermal spectrum from the NS as well as its overall cooling rate. In this paper, we 
consider the case when the rate limiting step for the H consumption is diffusion to the burning layer, rather than 
the local nuclear timescale. This is relevant for NSs with surface temperatures in excess of 10^ K, such as young 
(< 10^ yr) radio pulsars and accreting NSs in quiescence. When downward diffusion is the limiting rate in DNB, 
the rate of H consumption is suppressed by 1-2 orders of magnitude compared to a DNB estimate that assumes 
diffusive equilibrium. In order to apply our ongoing study to young neutron stars, we also include the important 
effects of strong magnetic fields (B ^ 10'^ G). In this initial study of magnetic modifications to DNB, we find that 
the H burning time is lengthened by 2-3 orders of magnitude for a 10'^ G field. However, even for NSs with dipole 
field strengths of 10'^ G, we find that all of the H can be burned before the pulsar reaches an age of ~ 10^ yr, thus 
potentially revealing the underlying proton-capturing elements. Finally, we conclude by providing an overview of 
what can be learned about fallback and pulsar winds from measuring the surface composition of a young NS. 

Subject headings: diffusion - magnetic fields - nuclear reactions, nucleosynthesis, abundances - pulsars; general 
- stars: abundances, interiors - stars: neutron 



1. Introduction 

The surface compositions of young NSs are now regularly 
probed by Chandra and XMM-Newton observations (see Pavlov, 
Zavlin & Sanwal 2002a for a review; also see Sanwal et al. 
2002; Mereghetti et al. 2002; Bignami et al. 2003). Since 
the amount of matter in the photosphere is small lO'^^M©), 
there is no reliable way to predict the composition of these outer 
layers from supernova theory. Fall-back during the supernovae 
event (or later) could provide a large range of possible elements 
(Woosley & Weaver 1995), ranging from pure H (due to spal- 
lation of heavier elements; Bildsten, Salpeter & Wasserman 
1992), mid- weight material (due to a "soft-landing" of outer 
shells falling back), or iron (if no fallback occurred). 

The current X-ray observations of young radio pulsars sug- 
gest that those younger than ~ lO^^yrs possess magnetic H 
or He atmospheres (Pavlov, Zavlin & Sanwal 2002a), whereas 
pulsars older than 10''~''yrs possess heavier element atmo- 
spheres (see Pavlov et al. 2002b and references therein; also 
see Kaminker, Yakovlev & Gnedin 2002). This suggests a pos- 
sible evolution of H or He to more blackbody-like elements on a 
timescale of 10'*~^yrs, and motivated our work (Chang & Bild- 
sten 2003; hereafter CB03) on a mechanism of nuclear evolu- 
tion called diffusive nuclear burning (DNB) (first mentioned by 
Chiu & Salpeter 1964 and later calculated by Rosen 1968) that 
provides for a surprisingly rapid depletion of H. 

The physics of DNB is very simple. At the photosphere, the 
local temperature (T^ ~ 10^ K) and density (p ^ Igcm"-') are 
far too low to drive nuclear burning. However, 1 meter under- 
neath the photosphere, the density and temperature are much 
higher [p ~ lO^'^gcm""*, T lO^K) so any surface H that dif- 



fuses down to this region will be rapidly captured onto a heavier 
nucleus. For example, for H capture onto Carbon, the '^C pro- 
duced in this reaction (following the f3 decay of '^N) sinks into 
the NS. This H depletion at depth then drives a H current from 
the surface that, over time, can completely deplete the H at the 
photosphere. The surprising realization is that the consumption 
of H is set by burning which occurs in the exponentially sup- 
pressed diffusive tail, and can easily consume all of the H on a 
NS in timescales of ~ 10^ years. 

In CB03, we found that at the characteristic depth, jbum (i-C- 
the burning layer), where the H is rapidly consumed, there was 
adequate time for diffusive equilibrium to be established for 
Tc < 5 X 10^ K for H on C. Hence at these low T^'s (and there- 
fore low Tg's), the rate limiting step is nuclear processes rather 
than diffusion. However, at higher T^'s typical of young NSs, 
the nuclear timescale, Tnuc, rapidly drops below the diffusion 
timescale, Tdiff, so that diffusive equilibrium no longer holds 
and downward diffusion to the burning layer is the rate limiting 
step. The simple estimate for this transition to diffusion-limited 
DNB is when Tnuc ^ Tdiff, where rjiff = li^/V is the time for H 
to diffuse an ion scale height, h. We find that incorporating the 
effects of diffusion into DNB decreases H consumption by 1-2 
orders or magnitude or greater at high Tc compared to a naive 
calculation assuming diffusive equilibrium. In addition, the de- 
pendence of the rate of H consumption with Tc changes when 
diffusion is the rate limiting step. Because of the size of the 
effects of diffusion at high 71 on DNB, we only use an approx- 
imate treatment of diffusion in our estimates of DNB rates. 

In addition, we also include the modification of the thermal 
profile due to the opacity changes induced by such high B fields 
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(Ventura & Potekhin 2001; Potekhin & Yakovlev 2001), since 
nearly all young neutron stars are also strongly magnetic {B ~ 
lO'^G). We find that strong magnetic fields alter the burning 
rate due to the modification of the thermal profile by 2-3 orders 
of magnitude. 

Applying our estimates to young NSs, we find that the age 
after which DNB is no longer active is roughly 10^~* years for 
standard cooling (and much sooner for rapid coohng). Prior 
to the cessession of DNB, H is always depleted on a timescale 
shorter than the present age. Following Michel (1975), we show 
that pulsar wind excavation of the NS surface is potentially a 
powerful mechanism for composition evolution. In particular, 
the pulsar wind can excavate to depths that exceed the maxi- 
mum thickness allowed for a H or He layer due to thermonu- 
clear constraints. 

This paper is structured as follows. In § 2 we discuss the ba- 
sic equations and their modification in the limit where diffusion 
is the rate limiting step. We discuss the effects of the magnetic 
field on the envelope's thermal structure and DNB in § 3, which 
then allows us to apply our work to young neutron stars in § 4. 
We also speculate on the possible excavation of material of NSs 
due to ion loss in a pulsar wind. We conclude in § 5 by sum- 
marizing our results and providing our views on what can be 
learned from measurements of young NS surface composition. 

2. Diffusion Limited DNB 

We only consider the top 10^ cm of the envelope, where the 
density is < lO'^'gcm"-'. The thickness of this layer is <S R, so 
we assume a plane parallel atmosphere with constant downward 
gravitational acceleration, g. Consider an envelope consisting 
of protons, electrons and ions that are coupled to each other and 
obey charge neutrality (jie = J2i^i'^i)- The equations of hydro- 
static balance are, 
dP 

= -np [mpg-eE+repneiVp-Ve) + ripniZfiVp-Vi)] ,(la) 

dPi r 2 

— = -Hi [Aimpg-ZieE + repHeZf (v, -Vg)-l- 

ripnpZf(vi-Vp)] , (lb) 

dP 

-^=-11, [m,g + eE + r,i,niZf(v,-Vi) + r,i,np(Ve-Vp)] ,(lc) 

where P,, n,. A,, Z,, v, = are the pressure, number den- 
sity, atomic number, charge and relative velocity of the f th ion 
species in Eulerian coordinates and E is the upward pointing 
electric field. General relativistic effects have been ignored. For 
simplicity, we presume the ions obey an ideal gas equation of 
state, neglecting effects such as plasma nonideality. We also ig- 
nore thermal diffusion effects, which are expected to play only 
a minor role in determining the structure of the NS envelope as 
the temperature scale height is larger than the ion scale height. 
In addition the thermal diffusion coefficient is small and hence 
thermal diffusion effects are expected to be minor compared 
to ordinary diffusion or gravitational settling in dense plasmas 
(Paquette et al. 1986 in the case of trace fuUy-ionized C in fuUy- 
ionized He). Plasma non-ideality effects are more significant 
especially for determining the structure of an envelope with two 
or more species with the same A/Z (i.e. He, C, N, O) and is an 
interesting question that needs to be addressed in the future. In 
the dense regions of the plasma, plasma non-ideality effects are 
expected to come in as Pcouiomh/P, where Pcouiomb would be the 
Coulomb pressure and P is the local pressure, which could be 
large (as large as 20%) in the region of interest. Plasma non- 
ideahty effects could be particularly significant in the case of 



highly magnetized plasmas (B > 10" G) as the Fermi pressure 
is decreased by the presence of a quantizing magnetic field (Lai 
2001). Hence our model is very simplified, but captures the 
essential physics of this problem. All vectorial quantities are 
positive radially. The coupling coefficient, rji, is due to scat- 
tering of species j by species 1. In CB03, we consider the case 
where these velocities are small, so that these drag terms could 
be neglected. However, this is not the case here. 

We consider the case where the protons are a drifting trace 
and the ions and electrons are a fixed background so Ve » v,- » 
and He « Zim » rip. Equations (1) simplify to 



dP, 



= -Hp {mpg-eE+ripriiZi Vp) , 



dr 

dPi 
'd? 
dP^ 

dr 



= -ni [AiMpg-ZieE) . 



= -«e {nieg + eE). 



(2a) 
(2b) 

(2c) 



We have dropped the drag terms in the ion and electron hydro- 
static balance equations since the drag force exerted on them by 

the protons is <C mpg. The drag term in the proton hydrostatic 
balance equation is only relevant in the presence of a nontrivial 



proton current, Jp 



The drag coefficient is simply re- 



lated to the diffusion coefficient, V, by taking the definition of 
the current, 

(3) 



Jp = -T^-^ +npWp = HpVp, 



where Wp is the proton drift velocity, and relating it to our equa- 
tions of hydrostatic balance. To begin let us presume that no 
current exists (/,, = 0) so that equation (3) becomes din Wp/dr = 
Wp/V. For an isothermal atmosphere the proton hydrostatic 
balance equation is 



dhiHp -nipg + eE 



(4) 



dr ksT 

which when equated to that from the diffusion equation gives 
Wp/'D = {—mpg + eE)/kBT. For Jp ^ 0, the proton equation of 
hydrostatic balance is 

dXnUp _ -mpg + eE-ripHiZfvp 



dr ksT 
whereas the diffusion equation becomes 
dlnnn w„-v„ 



(5) 



(6) 



dr V 

Taking our earlier expression for Wp and equating equations (5) 
and (6) yields the relation, 

riiZfV 



(7) 



between the drag and diffusion coefficients. 

In the limit of ideal gases, when Coulomb corrections to the 
ion equation of state are small (F = Z^e^ /{akBT) < 1, where a 
is the interionic spacing), we parameterized T) using the results 
of Chapman & Cowhng (1952) (also see Alcock & Illarionov 
1980). In the liquid regime (1 < F < 175), where most of the 
burning occurs, we begin with the fit of Hansen, McDonald & 
Pollock (1975) of the dimensionless self diffusion coefficient, 
V* = V,/{ujpa^)2.95T-''l'^, where LOp is the plasma frequency 
and a is the mean ion separation. We compute the interspecies 
diffusion coefficient from the self-diffusion coefficient in the 
same manner as Brown, Bildsten & Chang (2002). We assume 
the Stokes-Einstein relation, Vs/ksT = (47rai77)~\ where ai is 
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the charge neutral sphere of the background species and 77 is the 
viscosity. The diffusion coefficient for a trace in a background 
is VlksT = (47ra2??)~\ where a2 is the charge neutral sphere 
around the trace. Hence we find V = Vsai/a2, which gives 



V: 



^0.17-1.3 

yl.i'vO.i Jj.o ' 
z-j 



(8) 



where 1 and 2 are the background and trace ion respectively, 
Tfi = r/10''K and = p/ 10^ gem*-'. The burning layer is usu- 
ally in the liquid region unless the effective temperature is ex- 
tremely high or low. The self-diffusion coefficient fitted by 
Hansen et al. (1975) is off by 50% at F » 1, compared to the 
more complete calculation of Paquette et al. (1986), but is accu- 
rate to within a few percent for F > 10. The work of Paquette 
et al. (1986) presents a more consistent way of bridging the 
transition from the dilute approximate of Chapman & CowUng 
and the strong coupling results of Hansen et al. (1975), but for 
our estimates equation (8) is adequate in the liquid regime. At 
extremely low Tg's, the burning layer crystalhzes (F > 175 for 
a one-component plasma; Potekhin & Chabrier 2000) and halts 
diffusion. 

The thermal structure is determined by the heat diffusion 
equation for a constant flux, 

dT 

d7 



3«P „4 
16^3^"' 



(9) 



where k is the opacity. The radiative opacity is determined 
mostly by free-free absorption, where the Gaunt factor, ^ff, is 
fitted by Schatz et al. (1999). We use the conductive opacities 
given by Baiko et al. (1998) with analytic formulae given by 
Potekhin et al. (1999a). For NSs with low surface magnetic 
fields (B < 10' G), we use the Paczynski (1983) electron equa- 
tion of state (which is accurate to within 5% for ideal electron 
plasmas). 

The continuity equation for the protons is 



drip dJp 
dt dr 



' p.nuc 



(10) 



where Tp.nuc is the local proton lifetime against nuclear cap- 
ture. As we have argued in CB03, the timescale for chang- 
ing the local number density is the timescale for changing the 
Hydrogen column, which is much longer than the local nu- 
clear timescale. Hence, we take a steady state approximation 
{dnp/dt = 0), which provides the current induced by DNB. We 
combine equations (3), (5) and (10) to obtain 

Pi- T-r 

f^B^ 'p,nuc 



where Ip = ksT / (eE-nipg) is the proton scale height. 

Another way to derive equation (1 1) is to combine the proton 
continuity equation (eq. [10]) with the definition of the current 
(eq.[3]) to obtain 



d Hp Wp drip 
'~d^^^~d^^VT, 



= 0, 



(12) 



p,nuc 



where for the purposes of this discussion we presume V and 
Wp are constants. Our full analysis includes the position depen- 
dence of V and Wp in equation (12). However, we will argue 
that only two terms really matter in the final analysis. Taken 
with equation (7), we see that equation (12) is equivalent to 
equation (11). 

Care must be taken when solving equation (12) since rp_nuc 
has a powerful temperature dependence. We write rip = TpUpf) 



in terms of a scaling factor J^p and ripp, the solution in diffusive 
equiUbrium. After dividing by np^, equation (12) becomes, 



IlL (- 

"p.o V 



V 



'p,0 



dr^ dr 
d^J^, 



V- 



V 



dr 

dTp dlnnpfl _ 



(13) 



dr^ ^ dr dr 

The first term in parentheses is the homogeneous diffusion equa- 
tion and hence it is zero, while the second term, -dXnnpaj dr+ 
l~p , is the ion equation of hydrostatic balance and is solved by 
Wp^O) niaking this term zero as well. Thus the fuU diffusion equa- 
tion reduces to 



dr^ dr dr 



p ^dJ^p dXnripfi 



(14) 



p,nuc 



which we solve with the boundary conditions J^p^^l and dTp/dr f 
far above the burning layer. Since Tp,nuc is a nuclear timescale, 
which falls very rapidly as we move into the deeper, hotter re- 
gions of the envelope, the nuclear scale height, /^y^. = Prp nuc, 
rapidly decreases into the envelope. Hence the second term 
in equation (14) rapidly becomes insignificant as the proton 

scale height, Ip = [d In n^.o / dr^ becomes /nuc at the burning 
layer. Thus for all practical purposes, the second term could be 
dropped and we would have a diffusion equation for Tp with 
a nuclear driven source. Therefore, Tp defines the departure 
from diffusive equilibrium at a given depth as a result of nu- 
clear burning. If we include the position dependence of V and 
Wp, this would not change our result significantly since inclu- 
sion of those terms would introduce additional terms with Ip 

and the thermal scale height, Ij = (dlnT / dr^ ' , both of which 
are » Znuc at the burning layer. The total diffusive current is 
given by 



(,H = ■ 

Tcol 



^ p^pfl 



dr: 



. f ^pnpfi y 

J ^p.nuc P 



dlny, (15) 



where y = P/g is the total column density, p is the local density 
and Tcol is the H colunm lifetime. 

Figure 1 shows the thermal and compositional structure of 
the envelope (eq. [2a], [9] and [14]) and the total burning rate, 
(//. We also solve the model assuming full diffusive equilib- 
rium (CB03). The burning cuts off the H number concentra- 
tion at large columns. The exponential increase in burning 
rates modifies the solution for the total diffusive current by a 
factor of two compared to the case of full diffusive equilib- 
rium. We showed that the condition of DNB in diffusive equi- 
librium is satisfied when the temperature in the burning layer 
T(p) < 4.55 X lOV;" " K (CB03). The particular model that we 
have chosen (Te = 1 .5 x 10'' K) is in the diffusion limited regime 
and illustrates the reduction in the burning rate when diffusion 
becomes the rate limiting step. Compared to the calculation 
assuming diffusive equilibrium, incorporating diffusion effects 
decreases the burning rate, (h, by a factor of 25. This difference 
between a naive calculation which assumes diffusive equilib- 
rium with one that includes the effects of diffusion gets even 
larger at high T^'s. 

The effect of incorporating the effects of diffusion into the 
DNB calculation is illustrated by Figure 2, which shows the H 
column lifetime, Tcoi, as a function of core temperature, Tc. In 
nuclear limited DNB, the column burning rate follows the scal- 

ings derived in CB03 for H/C, TcoI oc y^^'^^T''^-^ (for a range 
of Tc between 2-4 x lO^K, but as noted in CB03, the scal- 
ing varies depending on the temperature we expand about; for 
lower temperatures, < 2 x 10^ K, the exponential scaling with 
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Tc is -12 to -13). As we transition to diffusion limited DNB, the 
scaling with central core temperature becomes more gentle as 
shown in Figure 2. The scaUng of Tcoi with total H column, 
jH, remains the same as in nuclear limited DNB, since the H 
number density in the burning layer scales like a y~j^' (CB03), 
where (5, = (Z, + l)/A,-2 « -1.42 for H on C, and is totally 
independent of the details of the burning. 

3. Effects of Strong Magnetic Fields on DNB 

In order to apply our work to radio pulsars, we must include 
the effects of strong magnetic fields (B ^ lO'^ G) on the micro- 
physics of the NS envelope. Magnetic fields of > 10" G affect 
the thermal structure of the envelope by modifying the electron 
equation of state and opacity (for a review see Lai 2001 and ref- 
erences therein). The changes to the electron equation of state 
modify the electric field, which modifies the diffusion problem. 
However, we neglect magnetic effects on the ions since mag- 
netic fields > 10'^ G would be required to modify the ion EOS. 

To model the full range of the electron EOS and opacity, we 
adopt the formalism of Potekhin & Yakovlev (2001) and refer- 
ences therein to relate the electron pressure and electron num- 
ber density to the chemical potential. We calculate the chem- 
ical potential from the local pressure by modifying the itera- 
tive procedure described by Potekhin & Yakovlev (1996). The 
free- free opacity in a magnetic field is from equation (21) and 
(22) of Potekhin & Yakovlev (2001). For electron conduction 
opacities, we use the results of Potekhin (1999). Starting from 
the photosphere at r = 2/3, our numerical integration of the 
constant flux equation (9) for the magnetized envelope agree 
to within 1% with the thermal profiles (i.e. T{p)) given by 
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Fig. I . — Total burning rate of H taking into account p-p capture and p + 
'^C capture. The bottom panel shows the number fraction, fy = nu/ntot, and 
temperature (dotted Une) as a function of column. This particular model has a 
total H colunm of 100 g cm~^ with an effective temperature of 1 .5 x 10* K. For 
the burning rate and number fraction, the dashed Unes are the solutions taking 
diffusion limiting into account. The solid lines are the solutions assuming 
diffusive equilibrium (CB03). Assuming diffusive equilibrium, the burning 
rate of this model is Of ^ 6300 gcm"^ yr"' . When diffusion Umiting is taken 
into account, the burning rate is Ch ^ 400gcm~^yr~^ 



Potekhin et. al. (2003) for accreted magnetized envelopes. 

The thermal structure of a Hydrogen on Carbon magnetized 
envelope is given in Figure 3 for magnetic fields of various 
strengths and orientation, which we define in terms of 9, the an- 
gle of the magnetic field relative to the radial direction. These 
NSs have surface gravities of ^ = 2.43 x lO^'^cms"^, for a 1.4 
Mq, 10 km NS with an effective temperatures of Tg = lO^K. 
We also show the nonmagnetic case (dotted Une) for purposes 
of comparison. Compared to the nonmagnetic case, the mag- 
netic envelopes are different in several respects. For a magnetic 
field with a significant radial component ( i.e. cos9 = 1 or 0.5), 
the base temperature is lower for a given effective temperature 
compared to the nonmagnetic case, while for more horizontally 
oriented magnetic fields, the base temperature is higher for a 
given effective temperature, highlighting the effect of magnetic 
field orientation on the opacity of the envelope. The transi- 
tion from primarily radiative to primarily conductive heat tran- 
port, after which the envelope becomes isothermal, takes place 
at larger column and hence higher density compared to the non- 
magnetic case. For the magnetic cases, orientation rather than 
field strength is a stronger determininant of the thermal profile, 
because the opacities along the field lines are a 4-7 orders of 
magnitude less than opacities across the field lines (for a re- 
view see Potekhin & Yakovlev 2001 and references therein). 
Therefore the temperature distribution is highly latitude depen- 
dent, with the star hotter at the magnetic poles and cooler at the 
magnetic equator. 

The electric field determines the Hydrogen concentration, 
fn, in the Carbon layer (CB03) and is a function of the adiabatic 
index, Xp = (dlnPg/dlnp) oc (dlnPe/dlnne). In the absence of 
magnetic effects, the adiabatic index undergoes a smooth tran- 
sition from = 1 in the ideal gas case to Xp = 5/3 in nonrela- 
tivistic degeneracy to = 4/3 in relativistic degeneracy. How- 
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Fig. 2. — Lifetime (solid hne) of Hydrogen with a fixed total column, y = 
lOOgcm"^ on Carbon and effective temperature (dashed line) as a function 
of core temperature. The model with effective temperature, Te = 10* K, or 
core temperature, Tc = 5.5 x 10^ K, represents the transition from the nuclear 
limited to diffusion limited regimes of DNB for H on C. 
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ever, in the presence of a quantizing magnetic field, the adia- 
batic index also undergoes de Haas-van Alphen oscillations in 
the weakly quantizing regime (Potekhin, Chabrier & Shibanov 
1999b; Lai 2001) that also appears in the electric field. 

Figure 4 shows the adiabatic index, Xp, and electric field as a 
function of density for the nonmagnetic case (fi = 0), = 10^^ G 
and B = 10^^ G. The de Haas-van Alphen oscillations are evi- 
dent as the electron EOS adiabatic index climbs from the value 
for the ideal gas equation of state to that of a 1-D Fermi gas, 
Pe cc p^. As the Fermi energy, Ef exceeds the integer values 
of the cyclotron energy, the equation of state softens, leading to 
alternating rising and falling of the adiabatic index and electric 
field. These oscillations are smeared out at the larger densities 
and temperatures deeper in the NS as Ef greatly exceeds the 
cyclotron energy and we recover the nonmagnetic solution. 

For simpUcity, we have ignore the effects of a strong mag- 
netic field on diffusion. However, we expect the effects of mag- 
netic fields on diffusion in the burning region is of little impor- 
tance at typical pulsar fields. The following estimate illustrates 
this point. The mean free path for proton diffusion across mag- 
netic fields is roughly given by the proton-gyro radius which is 
given by = y^WJJ^/ujB « 10~**(r/10'*K)i/2B-i cm, where 
ojb = eB /{nipc) is the cyclotron frequency. The mean free path 
for proton diffusion is given by the mean ion spacing a =[3p/ 



-1/3 



cm. Since the ion spacing 



is typically an order of magnitude smaller than the proton-gyro 
radius, we expect the effects of typical pulsar magnetic field on 
diffusion in the burning region to be small. Magnetic effects on 
diffusion is significant in the low density regions of the enve- 
lope (p ~ lOgcm"^), but changes to the diffusion coefficient in 
this region is of little importance to diffusion-limited DNB. For 
magnetar strengths, this is no longer the case and the effects of 
magnetic fields on diffusion must be taken into account. 
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Fig. 3. — Thermal structure of a magnetic neutron star envelope for various 
values of field strength and inclination where cos 6 = 1 is the magnetic pole. 
The effective temperature is fixed at r<, = 10* K and a H/C envelope is assumed 
with the H/C boimdary aiy = 10^ g cm"^. The thermal structure for a B=0 
envelope is shown by flie dotted line. 
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Fig. 4. — Electric field and adiabatic index Xp = dlnPe/dlnp for a pure 
H envelope for various B-fields. The shift of the degeneracy point for the 
higher B-field case is readily apparent for the B = lO'^G case. The initial 
approach of the adiabatic index to three for the higher B-field case and the 
response of the de Haas-van Alphen oscillations to B-field strength is also 
shown. As expected all three cases asymptotically approach the relativistic 
degenerate limit of eE/nipg = 1 andd\nPe/d\np = 4/3. 
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Fig. 5. — Total burning rate of H taking into account p-p capture and p 
+ '^C capture for a magnetized envelope with B = lO'^G at the magnetic 
pole (cos 6 = 1) and unmagnetized envelope with the same core temperature, 
Tc = 5.5 X 10^ K, and the same total H column, yn = 100 gcm"^. The bottom 
graph shows the number fraction and temperature as a function of column for 
the magnetized (solid and dotted lines) and unmagnetized (dashed and long- 
dashed lines) envelopes. The nonmagnetic envelope has a column lifetime of 
TH ~ 7yr, while the magnetic envelope has th ~ 200yr at B = lO'^ G. 
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In Figure 5, we plot the burning integrand, d(H/(dlny), for 
a model envelope with = 1.125 x 10^ K, B = lO^^G at the 
magnetic poles (cos^ = 1) and compare this against unmagne- 
tized model with the same core temperature and same total H 
column, yh- Compared to the unmagnetized model, the mag- 
netized model atB= 10'^ G burns at a rate that is an order of 
magnitude less compared to the unmagnetized model. The rea- 
son for this lowered burning rate is a combination of lowered H 
number density and lowered temperature at the burning layer. 
The bottom plot of Figure 5 illustrates both points. At the burn- 
ing layer (which is roughly the same for the magnetized and un- 
magnetized cases), the magnetized case for B = 10^^ G has both 
lowered temperature and lowered H number density. Another 
crucial difference is that, despite possessing the same core tem- 
perature, Tc = 5.5 X 10^ K, the magnetic envelope (B = 10'^ G) 
has a higher effective temperature (J,. « 1.1 x lO^K) than the 
nonmagnetic case (Tg w 10^ K) due to the lower radiative opac- 
ity of highly magnetized plasmas. 

Because the H number density or concentration and temper- 
ature profile are both modified by the presence of a quantizing 
B-field, the magnetic effects on DNB are nontrivial. The B- 
field always lowers the temperature at a fixed column for fixed 
Tc as shown in the bottom plot of Figure 5. The higher effec- 
tive temperature of the magnetized envelope is wholly due to 
the increase in depth where the photosphere sits. If the burning 
layer was at a fixed column, this would suggest that increasing 
the B-field tends to decrease the DNB burning rate. In Figure 
5, the difference in temperature at a given y between the mag- 
netized envelope and B = Q envelope is about 12%, which gives 
a factor of five difference in the capture rates between the two. 
Also the transition region between the radiative outer envelope 
and isothermal inner envelope and core (the sensitivity strip) 
moves toward greater column or density with increasing B-field 
(Ventura & Potekhin 2001). This would also decrease the DNB 
burning rate as the burning layer moves to deeper depth and 
hence decreasing H number density (at least in the case of nu- 
clear limited DNB). Also as illustrated in Figure 4, the electric 
field in the 10'^ G magnetized envelope is on average higher 
leading up to the burning region (p ~ 10"^^ gcm"^) than in the 
nonmagnetized envelope. This results in a lower number den- 
sity at a given column in the magnetized envelope (B= lO'^G) 
compared to the nonmagnetized envelope, as shown in the plot 
of number concentration as a function of colunm in Figure 5. 
Hence, we find for a B = 10^^ G field, the DNB burning rate is 
reduced by a combination of a modified thermal profile and a 
modified electric field, but for stronger fields the situation may 
be significantly different and is a topic of further study. 

We plot the H colunm lifetime as a function of the central 
core temperature for differing inclinations ofaB=10'^G with 
jH = lOOgcm"^ in Figure 6. For purposes of comparison, we 
also show the nonmagnetic case (thick solid line). The low and 
high temperature scalings of the magnetized envelope are sim- 
ilar to the scalings the nonmagnetized envelope, but the over- 
all normalization and transition temperature are modified by 
the magnetic field. The H column lifetime is greater for the 
magnetized envelopes at B= lO'^G than for the unmagnetized 
envelopes for at a given jh and Tc, but the magnitude of the 
change is highly dependent on the orientation of the magnetic 
field. In the bottom plot, we show the effective temperature as 
a function of Tc. The variance in column lifetime and effective 
temperature are mainly due to the modification of the temper- 
ature profile by the magnetic field especially in the conductive 



part of the envelope. 

For a dipole geometry, Figure 7 shows the latitude depen- 
dence of the effective temperature and column lifetime for fixed 
Tc, dipole field strength {Bp = 10'^ G) and total H colunm of 
yu = lOOgcm"^. The effective temperature and burning rate are 
greatest at the magnetic poles and smallest at the equator. At 
the magnetic equator the burning rate is reduced by about two 
orders of magnitude, while the effective temperature is reduced 
by about a factor of three. This change in effective temperature 
and column lifetime is mainly due to magnetic field inclination 
changing the region where the burning layer sits. From this 
plot, roughly 60% of the NS surface has the same burning rate 
(within 50%) and effective temperature (within 30%). Since 
H is preferentially depleted near the magnetic poles, this may 
induce a surface current of H from the equator to the poles, a 
problem that merits further study. 

Much caution must be taken when looking at these results 
seriously especially in the superstrong B-field case since the ef- 
fects of plasma non-ideahty become more and more pronounced. 
Changes in the ion and electron equation of state due to Coulomb 
interactions changes the electric field in these envelopes, possi- 
bly altering the diffusion problem drastically. 

4. Surface Composition Evolution for Young 
Radio Pulsars 

The clear application of our theoretical work is to young ra- 
dio pulsars, which can have both the nuclear physics of DNB 
active, but also, potentially mass loss via pulsars winds. Given 
the uncertainty in fallback models, the ease (in terms of avail- 
ability) with which H can be accreted or produced by fallback 
(i.e. spallation; Bildsten et al. 1992), and the tiny amount of 
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Fig. 6. — Hydrogen column lifetime and effective temperature for H on 
C envelopes as a function of core temperature Tc for different magnetic field 
inclinations, where cos 6 = 1 are the magnetic poles. The total H column is 
fixed at JH = 100 gcm~^. The magnetic field strength is B = lO'^G, typical 
of young radio pulsars. For a given core temperature the effective temperature 
of the magnetized object is greater for radial inclinations of the magnetic field, 
but the column hfetime of H is greater for the magnetized envelopes, regardless 
of inclination, than the nonmagnetic case (thick solid Une). 
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material (~ 10" M©) that comprises the NS photosphere, NS 
atmospheres should comprise mainly of Hydrogen. Further- 
more, the strong gravity ~ 2-3 x lO'^'cms"'^) on the surface 
of a NS implies that the sedimentation timescale is of order one 
second (Bildsten, Chang & Paerels 2003) and hence the photo- 
sphere should be pure H. Hence we would expect that all young 
NSs possess surface H at early times, but when the surface H 
begins to evolve depends on when the fallback episode ends, 
which is very uncertain. 

4.1. Surface Evolution via DNB 

In Figure 8 we show the photospheric H lifetime (tcoI at the 
photosphere) curves (i.e. implied Tc for a given lifetime) for 
B = and B= lO'^G and proton capturing elements C and O. 
With the view that H could either be initially present or added at 
late times, we must compare these burning rates to the amount 
of time the NS spends with this Tc. These times are available 
from NS cooling calculations, and in Figure 8 we overlay a few 
representative cooUng tracks (long-dashed Unes). We include a 
standard cooling track which presumes modified URCA cool- 
ing on a 1.3 Mq NS and another with core proton superfluidity 
(Potekhin et al. 2003). D. G. Yakovlev kindly calculated a 
model for us with triplet-state neutron superfluidity in the core 
with a critical temperature of 8 x 10^ K. 

This figure makes clear that for H overlaying C, the burning 
time is always much less than the pulsar's age for the first 10^"^ 
years. After this active DNB phase, any new H could stably 
exist on the photosphere. For H on O, we find that the active 
DNB phase depends more sensitively on the NS cooling and 
varies from 10"*"^ yrs. This implies that H deposited on the sur- 
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Fig. 7. — Hydrogen column lifetime for yn = lOOgcm ^ and effective tem- 
perature for H on C envelopes as a function of inclination to the magnetic pole 

9, presuming a dipole field configuration, B = Bp cos 6 + sinO/lj, where 

Bp = 10^2 G, for three different core temperatures, Tc = 6.3, 8, 10 x lO'' K. The 
lifetime Tcoi changes by almost two orders of magnitude with latitude while 
the effective temperature changes only by a factor of a few. The variation in 
burning rate with latitude at fixed Tc is due to the changing temperature of the 
burning layer due to magnetic field inclination. 



face of a rapidly cooling NS after 10^ years could exist stably if 
the underlying material is Oxygen or heavier. No matter what, 
the active DNB phase does not naturally provide an explanation 
for the 10^^ yr timescale for surface evolution implied by the 

observations. 

The recent observations of mid-Z atomic fines on IE 1207.4-5209 
implies that any surface Hydrogen must have been depleted on 
at least a 10^ yrs timescale (Sanwal et al. 2002; Mereghetti et 
al. 2002; Hailey & Mori 2002; Mori & Hailey 2003). This 
pulsar is associated with the SNR, G296.5H-10; the age by as- 
sociation is « 7kyrs (Roger et al. 1988). The spectral features 
have been interpreted as either Oxygen Unes (B12 = 0.55-0.75, 
z = 0.06-0.21) or Neon lines (B12 = 0.8- 1.1, z = 0.62-0.86) 
(Mori & Hailey 2003). They have also been interpreted as He 
lines in a B ^ 10^^ G field, but this would require a special mag- 
netic geometry in order to agree with the spindown B-field of 
Bn « 2-4 (Pavlov et al. 2002b; Sanwal et al. 2002).^ 

Let us take the interpretation of these spectral features as 
Oxygen lines. Mori & Hailey's (2003) interpretation of the 
spectral features as Oxygen lines yielded B = 5.5-7.5 x 10^' G 
with a redshift of z = 0.06-0.21 . Though it is not self-consistent, 
we use the « 1 . 1 X 1 0^ K from the H atmosphere fit of Mereghetti 
et al. (2002). For Mori & Hailey's combination of redshift and 
B-field, we find that the H Ufetime at the photosphere varies 
from 2-20 kyrs for an underlying Oxygen layer. Hence, it is 
plausible that DNB has occurred on the surface of this NS to 
deplete any surface H (;/ initially present) to reveal the under- 




FlG. 8. — H column lifetime curves for H on C and H on O envelopes 
at the magnetic poles for different Bp = (thick solid line), 10'^ G (dotted 
line). Overlayed are cooling tracks (long-dashed lines) for various NS mod- 
els (Potekhin et al. 2003) including standard cooling (no superfluidity, mod- 
ified Urea cooling represented by 1.3M0 noSF), proton superfluidity (pSF) 
and triplet-state core neutron superfluidity (nSF) with critical temperature 
» 8 X 10* K. 

' More recent observations indicate that there are additional fines at 2.1 and 
2.8 keV (Bignami et al. 2003), which would render the intepretation of these 
lines as electron cyclotron lines, but these results require further study. The 
electron cyclotron lines would indicate a magnetic field sttength of B = 6 x 
10i'*(l+z)G. 
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lying proton-capturing element, Oxygen. 

4.2. Surface Composition Evolution via Pulsar Wind Mass 

Loss 

In addition to DNB, the surface evolution of magnetized ro- 
tating NSs may also be driven via outgoing pulsar winds. Though 
there is debate about whether nuclei are present in the pulsar 
wind (Goldreich & Julian 1969; i.e. see Ruderman and Suther- 
land 1975), it is still helpful to estimate how much material 
could possibly be excavated from an active pulsar's surface. 
Gallant and Arons (1994) analysis of the Crab pulsar winds 
showed that the ions leave that pulsar roughly at the rate of the 
Goldreich- JuUan current flow from the polar cap region, 

20^/icosj 



ec 



(16) 



where (i is the NS magnetic moment. We integrate this over 
time to find the total number of excavated ions, = J Nidt 
(Michel 1975). For a constant magnetic moment (i.e. no field 
decay), the time integral can be turned into an integral over 
the spin frequency using the magnetic dipole spin-down law, 
fj = -2/x^ri^sin^(?/3/c-', where / is the NS moment of inertia 
and 9 is the angle between the spin axis and magnetic moment. 
This gives Michel's (1975) estimate of the total number of ions. 



3/c^cos; / rii 

— —M^ 



(17) 



that have been excavated from the NS surface as it spins down 
from to f2 /. 

Though this has never been shown, we will presume that the 
long excavation timescale (of order the spin-down time) allows 
the polar "hole" left by the exiting ions to be filled in by flow of 
material from other parts of the NS. Taking the moment of iner- 
tia as, Ik QAME} (Ravenhall and Pethick 1994), the excavated 
column depth (presuming pure Hydrogen) is 



where we dropped the angular factors. Putting in fiducial num- 
bers, we find a total excavated column density 

10^0 G cm^ 



yexc ~ lO'g cm ^ 



M 



(19) 



,1.4M0/ V M 

which is comparable to the maximum amount of Hydrogen (Bild- 
sten & Gumming 1998) or Hehum (Gumming & Bildsten 2000) 
that can be present on young NSs with hot cores. 

Though the physics of the pulsar wind is by no means well 
understood, excavation has testable observational consequences 
for magnetized spinning NSs. Since >'exc is very sensitive to 
the NS's magnetic dipole moment, the possible repercussions 
are diverse. For conventional radio pulsars, the underlying ma- 
terial may be unveiled during the initial cooling phase (< 10*' 
yrs). For magnetars and AXP's, a large layer of overlying mate- 
rial can survive on the NS. However, for millisecond radio pul- 
sars, the amount excavated is so large that the rp-process ashes 
from the mixed Hydrogen-Helium burning during the accretion 
phase (see Strohmayer & Bildsten 2003 for an overview) might 
well be revealed. 

5. Conclusions 

Supernova theory cannot predict the surface composition of 
a young neutron star. The range of possible initial surface 



compositions is large, depends on fall-back physics and resid- 
ual burning at early times (i.e. < 1 - 10 years). Hence, we 
focused here on tracing the abundance evolution for different 
initial states in the hopes of connecting them to possible final 
states. This has allowed us to construct a guide to the observer 
(Table 1) on the astrophysical implications of secure detections 
of specific elements. Before explaining this observational guide 
in detail, we summarize our main physics results. 

We have extended our original study of DNB (CB03) to in- 
clude burning when diffusion is the rate limiting step and the 
important effects of surface magnetic fields in the context of a 
highly approximate model of the envelope and diffusion physics. 
When diffusion is the rate limiting step, we modified the ba- 
sic equations of the envelope by including an extra "diffusion" 
equation with an nuclear driven sink. By including magnetic 
effects, we can now apply our work to young radio pulsars. 
We have shown that H can be depleted in 10* years from the 
surface of a magnetized neutron star by diffusively burning in 
deeper C layers. The timescale of this burning depends pri- 
marily on three NS parameters: the core temperature, magnetic 
field and nuclear composition of the layers where the proton 
captures occur. For magnetic fields typical of young radio pul- 
sars (B ^ 10'"G) with an underlying C layer, we have shown 
that the active DNB phase is relatively long (up to 10* yrs). 
The active DNB phase gets shorter for heavier proton captur- 
ing elements. We appUed our work to the recent observations 
of magnetized Oxygen lines from IE 1207.4-5209 and showed 
that (if initially present) Hydrogen would be depleted from the 
NS surface on a timescale comparable to its age (< 7 kyrs), 
exposing a pure Oxygen layer. 

In section 4, we pointed out that the loss of ions in a pulsar 
wind may be a relevant mechanism for the surface evolution of 
magnetized, rotating NSs. Our estimate for the amount of ma- 
terial lost via the wind is comparable to the maximum amount 
of Hydrogen-Helium on the surface of the NS at birth. Though 
the mechanism for the wind is highly uncertain, this amount 
of excavated material may have significant observational con- 
sequences for young radio pulsars and MSPs. 

There is one remaining compUcation that needs to be re- 
solved, which is the impact of an intervening Helium layer on 
the H depletion rate. The He layer would need to displace the 
burning layer in order to significantly impact the rates. How- 
ever, the degeneracy of such a thick He layer allows the He to 
readily diffuse downward into the underlying same A /Z matter. 
Hence, both Coulomb and thermal effects on the ion equation 
of state must be included to resolve the miscibility of the He 
with the underlying layers. This problem is also apphcable to 
the study of stratification in cooUng white dwarfs and will be 
addressed separately. 

Given the sensitivity of DNB and excavation to the initial 
composition of the envelope and physics of the pulsar wind 
respectively, we can learn a tremendous amount from preci- 
sion surface composition measurements of young B ~ 10^^ G 
neutron stars. We have contructed a guide to understand such 
measurements in Table 1 . Given the efficiency of DNB at the 
high Tc typical of young NSs, detection of H on a young pulsar 
would indicate that either continual accretion of H rich material 
is occuring, or a thick He buffer that blocks DNB is present, or 
there are no proton capturing elements present underneath the 
H. If H were detected on an older pulsar (> 10* yrs), this would 
indicate that ion loss via pulsar wind excavation is minimal. 
Detection of He on a young radio pulsar (< 10* yrs) would in- 
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dicate that the initial H was depleted (via either DNB or excava- 
tion) or that the material which fell back on the NS shortly after 
birth made no H. For NS older than 10^ yrs (where here we 
presume this is roughly the pulsar spin-down time), detection 
of He would indicate that pulsar wind excavation is ineffective 
at removing He from the surface. Detection of CNO or heav- 
ier elements may indicate that the initial H/He (and CNO for 
heavy opaque material) pile was depleted via excavation (possi- 
bly DNB for H). Detection of CNO elements may also indicate 
that the outer shells had a "soft landing" fallback onto the NS 
after birth i.e. no spallation, if excavation is excluded. Finally, 
heavy opaque materials may also indicate lack of any fallback, 
if there is no excavation. 



We thank A. Y. Potekhin for clarifying discussions on the 
nature of the electron EOS and opacities at high densities and 
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Table 1 

Implications of Obslrvkd Slrfacl CoMPOsnTONS of Yolng 10'- G Nll-i ron Stars 



Measured 
Surface Composition 


AgeofNS 


Astrophysical Implications 


H 


< 10'' yrs 


continual accretion 

thick He buffer 

- or - 

no underlying p-capturing material 
same as for < 10^ yrs 

- and possibly - 

no pxcavation from niilsar wind 








> 10^ yrs 


He 


< lO*" yrs 


initial H depleted via DNB/excavation 

- or - 

snallalinn of fallbaclc matprial witboiit H nrodiiction 

same as for <C 10^ yrs 

- and possibly - 

no He excavation from pulsar wind 




> 10^ yrs 


C, N,0 




initial H/He depleted via DNB/excavation 
- or - 

"soft landing" of outer shell material (no spallation) 






Heavy opaque elements 




excavation of initial H/He/CNO pile 



- or - 

no SN fallback (Mfauback < IO-^'Mq) 



